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ABSTRACT 

Aims. IKT 16 is an X-ray and radio-faint supernova remnant (SNR) in the Small Magellanic Cloud (SMC). A previous X-ray study 
of this SNR found a hard X-ray source near its centre. Using all available archival and proprietary XMM-Newton data, alongside new 
multi-frequency radio-continuum surveys and optical observations at Ha and forbidden [SII] and [OIII] lines, we aim to constrain the 
properties of the SNR and discover the nature of the hard source within. 

Methods. We combine XMM-Newton datasets to produce the highest quality X-ray image of IKT 16 to date. We use this, in combina- 
tion with radio and optical images, to conduct a multi-wavelength morphological analysis of the remnant. We extract separate spectra 
from the SNR and the bright source near its centre, and conduct spectral fitting of both regions. 

Results. We find IKT 16 to have a radius of 37 ± 3 pc, with the bright source located 8 ± 2 pc from the centre. This is the largest 
known SNR in the SMC. The large size of the remnant suggests it is likely in the Sedov-adiabatic phase of evolution. Using a Sedov 
model to fit the SNR spectrum, we find an electron temperature kT of 1 .03 ± 0. 12 keV and an age of « 14700 yr. The absorption found 
requires the remnant to be located deep within the SMC. The bright source is fit with a power law with index T = 1.58 ± 0.07, and is 
associated with diffuse radio emission extending towards the centre of the SNR. We argue that this source is likely to be the neutron 
star remnant of the supernova explosion, and infer its transverse kick velocity to be 580 ± 100 km s -1 . The X-ray and radio properties 
of this source strongly favour a pulsar wind nebula (PWN) origin. 

Key words, galaxies: Magellanic Clouds - ISM: supernova remnants 
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emission is detected and much of the broad-band information is 
lost. In addition, the distance measurements to sources within 
the Galaxy are often uncertain by a factor of two, leading to an 



(120041) (hereafter VDH04). iFilipovic et al] (120081) presented 
XMM-Newton results from another three known and one can- 
didate remnant. These studies found that the majority of SMC 
* Based on observations obtained with XMM-Newton, an ESA sci- SNRs are hkely in the Sedov-adiabatic phase of evolution, 
ence mission with instruments and contributions directly funded by IKT 16 was discovered as a candidate supernova remnant in 
ESA Member States and NASA. the SMC bv llnoue et al.l d 19831) using Einstein, and its identity as 
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Table 1. Details of the XMM-Newton observations covering IKT 16 used in this study. 



Observation ID 


Start Date 


Filter" 


Pointing 


co-ordinates 


Useful exposure (ks) 


Off-axis angle 




(yyyy-mm-dd) 


pn/MOSl/MOS2 


RA (12000) 


Dec (J2000) 


pn 


MOS 1 


MOS2 


(arcmin) 


0018540101 


2001-11-20 


M/M/M 


00''59 m 26.8' 


-72° 09' 55" 


23. 4* 


18.1* 


18.0* 


9.3 


r\(\Q a o r\r\ 1 a 1 

0084200101 


TAAO A1 TA 

2002-03-30 


1/M/M 


00 56 41. / 


-no ^ni 11/' 

— 12 20 11 


8.8 


10.0 


10.3 


8.0 


0110000201 


2000-10-17 


M/M/M 


00*59 m 26.0* 


-72° 10' 11" 


14.3* 


16.7* 


16.8* 


9.1 


0212282601 


2005-03-27 


C/M/M 


00''59'"26.8' ! 


-72° 09' 54" 


0.0 


3.8 


3.8 


9.3 


0304250401 


2005-11-27 


M/M/M 


00''59'"26.8 ! 


-72° 09' 54" 


15.9 


17.4* 


17.4* 


9.3 


0304250501 


2005-11-29 


M/M/M 


00''59"'26.8' ! 


-72° 09' 54" 


14.9 


16.5* 


16.6* 


9.3 


0304250601 


2005-12-11 


M/M/M 


00''59'"26.8* 


-72° 09' 54" 


10.6 


16.4* 


16.2* 


9.3 


0500980201 


2007-06-06 


T/M/M 


01'' 00'" 00.0' 


-72° 27' 00" 


14.8 


23.2 


23.9* 


11.8 


0601210801 


2009-10-09 


T/M/M 


00''56"'15.5' 


-72° 21' 55" 


23.0* 


24.6 


24.6* 


10.4 



" - T = thin filter, M = Medium filter, C = Filter-wheel Closed. 
* - Observation used for spectral analysis. 



a s hell-type SNR was confi rmed through radio and Ha analy- 
sis dMathewson et alj|1984l) . VDH04 included this remnant as 
part of their study of SMC SNRs. Using the earliest XMM- 
Newton observation, they found it to be X-ray faint, and its prop- 
erties were hence poorly constrained. A region of harder X-ray 
emission close to the centre of the remnant, whose source was 
unclear, was also reported. 

Nine XMM-Newton observations covering IKT 16 now ex- 
ist, giving an eightfold increase in exposure time. In this paper 
we present a study of IKT 16 and the object located towards its 
centre using X-ray data from XMM-Newton, radio data from the 
Australia Telescope Compact Array (ATCA) and optical images 
from the Magellanic Clouds Emission Line Survey (MCELS). In 
$2l we outline the procedures used to reduce and analyse XMM- 
Newton data. We also discuss complementary radio and optical 
data. In §0 we present the results of spatial and spectral analysis 
of IKT 16 and the X-ray hard object inside it. In $H we discuss 
the derived properties of IKT 16 and speculate on the nature of 
the hard source. We present our conclusions in §0 

2. Observations and data reduction 

2.1. X-ray observations 

XMM-Newton has observed IKT 16 serendipitously with EPIC 
in full-frame mode on 9 occasions between 2000 and 2009, at 
off-axis angles between 8' and 12'. Eight of the nine observa- 
tions were used for analysis, with one (obsid 0212282601) dis- 
carded due to persistent high background. Details of the obser- 
vations are shown in Table [TJ 

Data reduction was based on SAS vlO.O.ffl Datasets were 
screened for periods of high background through the creation 
of full-field single event (PATTERN=0) lightcurves above 10 
keV (with an upper limit of 12 keV for pn data). MOS data 
were excluded when the count-rate in a 25 second bin exceeded 
0.4 count s , with pn data excluded above 0.5 count s _1 . These 
thresholds were selected based on quiescent levels in the obser- 
vations. After filtering, useful exposure times range between 9 
and 25 ks. 

For each observation, images and exposure maps were cre- 
ated in several energy bands: 0.3-0.8 keV, 0.8-1.2 keV and 1.2- 
2.0 keV for production of 3-colour X-ray images, and 0.3-1.0 
keV for a broader-band soft X-ray image. We used data from 
single and double pixel events in the pn camera (patterns 0-4) 
and single to quadruple pixel events in MOS cameras (patterns 
0-12). The selection FLAG=0 was used to exclude spurious data 



from bad pixels, hot columns and chip gaps, and the pixel size 
was set at 4" x4". A constant particle rate, estimated from the 
corners of each image not exposed to the sky, was subtracted 
from each image. The images for each camera were then co- 
added using the SAS task emosaic in order to account for dif- 
ferences in pointing direction between the observations, and the 
exposure maps were added in the same way. Finally, the mosaic 
image for each camera was divided by the relevant exposure map 
and the images from the three cameras summed to produce a flat- 
field image suitable for further spatial analysis. 

Spectral analysis was performed on 15 exposures (3 pn, 5 
MOS 1 and 7 MOS2), selected on the basis of exposure time and 
position on the detector. When studying extended sources, it is 
necessary to be careful in defining source and background ar- 
eas to extract. In this case there are two areas of interest: the 
central source (detected as a point source in our X-ray obser- 
vations) located at RA(J2000)=00 /l 58"T6.7 5 Dec=-72° 18' 06" 
and the remainder of the supernova remnant. The shock bound- 
ary of the SNR is well-described by a circular region centred on 
RA=00 /, 58 m 22.4 i Dec=-72° 17' 52" with a radius of 128 ± 8", 
with the error due to uncertainty in the best fit circle radius . We 
therefore defined two source regions for each observation: a 20" 
radius area centred on the bright source and a 128" radius area 
centred on our best estimate of the centre of the SNR, exclud- 
ing the previously defined source area. These extraction regions 
are shown on Fig.[TJ(top left). For the first region, a background 
region was extracted from a nearby position on the same CCD 
outside the SNR. For the second, a nearby background region 
of similar size to the source region was found such that the av- 
erage exposure time for source and background regions was the 
same. This was done to ensure that the contributions from the 
cosmic X-ray background (vignetted by the telescope) and par- 
ticle background (not vignetted to first order) remained in the 
correct proportion when applied to the source spectrum. The 
SAS tasks arfgen and rmfgen were used to generate appropri- 
ate Auxiliary Response Files (ARF) and Response Matrix Files 
(RMF) for each observation, and the counts in adjacent spectral 
bins were summed to a minimum of 20 counts per bin using the 
Ftoo^ grppha. 

2.2. Optical observations 

The Magellanic Cloud Emission Line Survey (MCELS) was 
carried out using the University of Michigan/CTIO Curtis 
Schmidt telescope. Using narrow band filters corresponding to 
rfo(6563A), [S n](6724A) and [Om](5007A) line emission, plus 



1 Science Analysis Software, http://xmm.esac.esa.int/sas/ 



2 http://heasarc.gsfc.nasa.gov/ftools/ 
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Fig. 1. Top left: Lightly smoothed RGB XMM-Newton EPIC image of IKT 16 (R=0.3-0.8 keV, G=0.8-1.2 keV, B=1.2-2.0 keV). 
The two regions marked are those used for spectral analysis for the SNR (white circle) and the bright source (black circle). The 
centre of the remnant is marked with a cross. The image is binned in 4" x4" pixels, and is displayed logarithmically with a dynamic 
range of 20. Top right: MCELS [S ii]/Hq; ratio of IKT 16. The image is binned in 2" x 2" pixels. The white circle indicates the 
extent of the SNR. Bottom left: Lightly smoothed soft X-ray image (0.3-1 keV) overlaid with high-resolution 20 cm (1400 MHz) 
ATCA radio contours. The contours are at levels of 0.07, 0.14 and 0.28 mJy beam -1 . The X-ray image is displayed logarithmically 
with a dynamic range of 20, and the X-ray colour bar is shown in units of log(pn+MOS) ct/pixel/20 ks. Bottom right: Soft X-ray 
map overlaid with 3 cm (8640 MHz) ATCA radio contours. The contours are at levels of 0.25, 0.5 and 1.0 mJy beam -1 . The beam 
sizes for the radio observations are shown in Table 



red and green continuum bands (used to subtract most of the 
stars from the images), images of diffuse structure across both 
Magellanic Clouds have been produced with an angular resolu- 
tion of 2.4" /pixel. An image of the [S n]/Ha ratio across the 
remnant is presented in Fig. Q] (t op right ) . Furt her information 
about the MCELS is given in lSmith et all d2000h . 

Integral Field Spectroscopy (IFS) was performed on 
3 November, 2010 at Siding Springs Observatory using the 
2.3-m Advanced Technology Telescope and its Wide Field 
Spectrograph (WiFeS). WiFeS provides a 25" x 38" field with 
0.5" per pixel spatial sampling along each of 25" x 1" slits. The 



output format match two 4096x4096 pixel CCD detectors in 
each camera individually optimized for the blue and red ends of 
the spectrum. The 900 second single exposure was made in the 
central region of IKT 16 at position angle (east of north) 90 de- 
grees under clear skies with seeingestimated at 1 ". The exposure 
was made in classical equal modqj using the RT560 beam split- 
ter and 3000 Volume Phased Holographic (VPH) gratings. For 



3 This is a principal WiFeS data accumulation mode in which data 
were accumulated in the red and blue cameras on the source for equal 
times. 
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these gratings, the blue (708 lines mm ') range includes 3200- d r 
5900A while the red (398 lines mrrT 1 ) range includes 5300- 
9800A. 

The data were reduced using the WiFeS data reduc- 
tion pipeline based on NOAO (National Optical Astronomy 
Observatory) IRAF software. This data reductio n package was 
devel oped from the Gemini IRAF package (McGregor et al.l 
2003). The pipeline consists of four primary tasks: wifes to 
set environment parameters, wftable to convert single exten- ^ 
sion FITS file formats to Multi-Extension FiTS ones and create g ^ 
file lists used by subsequent steps, wfcal to process calibration | § 7 
frames including bias, flat-field, arc and wire; and wfreduce to ° 
apply calibration files and create data cubes for analysis. 



2.3. Radio observations 

We have extracted all archival radio-continuum observations of 
IKT 16 using the Australia Telescope Compact Array (ATCA) 
comprising observation s at 4 radio wavelengths: 20, 13, 6 
and 3 cm for ATCA dPickel et alJl200il IFilipovic et all l2002t 
iPavne etail l2004t IFilipovic et alj 1 19971) . In addition we make 
use of the Molonglo Observatory Synthesis Telescop e (MOST) 
SMC survey at 36 cm (843 MHz) dTurtle et al.lll998l) . The beam 
size for these observations ranges from 7" to 45". IKT 16 
was also observed with the ATCA as a part of projects C281 
and C634. More information about observing procedure and 
other sources obser v ed during these sessi ons c an be found in 
Boiic ic etalJ d2007h . iBoiicic et all (l2010h and ICrawford et all 
(120081) . Soft X-ray images of IKT 16 overlaid with 3 cm and 
high-resolution 20 cm radio contours are shown in Fig. [1] and 
the integrated flux densities detected from our radio observa- 
tions of the SNR are shown in Table [2] The error in the mea- 
sured integrated flux density (for the whole system and the bright 
source) is estimated as a quadrature sum from the local noise 
level (0.1 mJy beam -1 ) and the uncertainty in the gain calibra- 
tion (10%). For the SNR region (excluding the bright source) the 
error is calculated by adding the errors for the other two extrac- 
tion regions in quadrature. 



3. Results 

3.1. X-ray 

3.1.1. Morphology 

Fig. □ shows a 3-colour X-ray (0.3-0.8 keV, 0.8-1.2 keV and 
1.2-2.0 keV) image of IKT 16, along with radio images from 
ATCA and optical data from MCELS. The X-ray images are 
lightly smoothed with a Gaussian kernel of width 4". We de- 
fine the boundary of the SNR to be where the surface brightness 
of the X-ray emission abruptly falls to the background level. 
Using this definition, the radius of the remnant is estimated to 
be 128" + 8", with the error resulting from a 2-pixel uncertainty 
in fitting a circular region to the edge of the X-ray emission. At 
the assumed distance of the SMC (60 kpc), this corresponds to 
a radius of 37±3 pc, consiste nt with bein g an older remnant in 
the Sedov phase of evolution dCoxll 19721) . The measured radius 
makes it the largest known SNR in the SMC. The shock bound- 
ary is well-defined in the north and west of the SNR, where the 
diffuse X-ray emission is strongest. We find small variations in 
the temperature measured across the remnant, with the region of 
strongest X-ray emission in the north slightly hotter than its sur- 
roundings. More tenuous emission is observed to the east of the 
remnant, and there is some evidence for greater extension of the 
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Fig. 2. Radio spectra of the entire SNR (blue), the bright source 
near the centre of the remnant (red) and the remainder of the 
remnant (green), constructed from radio observations at 5 fre- 
quencies using MOST and ATCA. The best fit power-law fits 
are shown, and the data is also shown in Table [2] 



SNR in this direction. Depending on the direction of measure- 
ment, the linear radius of the SNR varies between 122" (in the 
N-S direction) and 144" (NE-SW), implying that the remnant is 
not perfectly circular on the sky. These deviations are relatively 
small and our simplification of the geometry of the system does 
not have a large effect on the results we obtain. These variations 
do however lead to uncertainty in the position of the SNR cen- 
tre, which limits the precision with which we can measure the 
dynamics of the system. 

The bright source observed inside the SNR shell is sig- 
nificantly harder than the emission from the rest of the rem- 
nant, as seen in the 3-colour image and its X-ray spectrum (in 
Section 13.1.2b . We find no evidence for this source being ex- 
tended in XMM -Newton measurements. We therefore limit its 
maximum X-ray extent to the FWHM of the detector (6"), which 
is 1.7 pc at the distance of the SMC. Assuming that the SNR 
shell can be described as spherically symmetric, the position 
of the bright source is offset from the centre of the shell by 
30" ± 8" (with the error mostly due to uncertainty in the po- 
sition of the SNR centre). Not taking into account possible pro- 
jection effects, this corresponds to a linear distance of 8 + 2 pc. 
We note that the bright source has been pre viously identified i n 
a 9 ks Chandra observation (OBSID 2948: lEvans etalJl2010h . 
The Chandra X-ray position is consistent with that found here, 
but due to the extremely short observation time and large off-axis 
angle in that observation, the errors on the XMM-Newton posi- 
tion are smaller. The observation also suggests that the source 
may be slightly extended (with an extent of 1.1 "±0.6"). A longer 
on-axis observation is necessary to constrain this value precisely. 
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V 


A 


Beam Size 


S Total 


S Source 


SsNR 


Telescope/Project 


(MHz) 


(cm) 


(arc sec) 


(Jy) 


(Jy) 


(Jy) 




843 


36 


44.9x43 


0.087 


0.045 


0.042 


MOST 


1371 


20 


7.05x6.63 


0.071 


0.039 


0.032 


C281 


2378 


13 


7.02x6.58 


0.058 


0.038 


0.020 


C281 


4798 


6 


30x30 


0.041 


0.033 


0.008 


C634, Parkes 


8640 


3 


20x20 


0.030 


0.028 


0.002 


C634, Parkes 



3.1.2. Spectral analysis 

Spect ral analysis was performed using XSPEC v!2.5.1 dArnaudl 



1996). For each of the relevant observations, source and back- 
ground spectra were extracted as described in the previous sec- 
tion. All 15 spectra extracted for each of the two extraction 
regions were fitted simultaneously, allowing only for a renor- 
malisation factor between observations from different cameras. 
To account for photoelectric absorption by interstellar gas, two 
absorption components were used in fitting: a phabs compo- 
nent fixed at the foregroun d Galactic value of 6 x 10 20 erg s -1 
(Dickev & Lockman 1990) assuming elemental abundances of 
Wil ms et al.l (feOOOn Tand a variable absorption {vphabs) compo- 
nent to account for absorption inside the SMC. This second com- 
ponent ha s metal abundances fixe d at 0.2 solar, as is typical in 
the SMC (lRussell&Dopitalll992h . 

X-ray emission from the central region was fit with an ab- 
sorbed power-law, and emissi on from the remainder of the SNR 
was fit with a Sedov model fBorkows ki et al.ll200ll) . Due to the 
relatively broad point spread function of the XMM-Newton EPIC 
detector, we estimate that for a point source, ~25% of the source 
counts will fall outside an extraction region of 20" . Therefore, 
an extra power-law component is introduced to the SNR fit, with 
index fixed to that obtained through fitting the central region and 
normalization set to account for this spillover. Similarly, as the 
central extraction region encompasses 2.5% of the overall SNR 
region, an appropriately scaled Sedov component is added to the 
central source model, with parameters fixed to those derived for 
the SNR. 

The free parameters of the Sedov model are the mean shock 
temperature, electron temperature just behind the shock front, 
metal abundances, ionization age (electron density behind the 
shock front multiplied by the remnant age) and normalization. 
Providing that the SNR has spherical symmetry and is in the 
Sedov phase of evolution, the fits obtained can be used to esti- 
mate several physical parameters. Using the distance to the SNR 
(D) and radius (R in m, assuming a distance of 60 kpc), vol- 
ume of X-ray emitting material ( V in m 3 ), normalization derived 
from the XSPEC fit (AO, shock temperature (T$ in keV) and 
baryon number per hydrogen atom (r m w 1.4, assuming a he- 
lium/hydrogen ratio of 0.1), it is possible to derive the electron 
density (n e ), age of the remnant (td yn ), total emitting mass (M), 
initial explosion energy (Eq) and ionization age (I t ) through the 
follow ing system of equations (VDH04, from Borkowsk i et al.l 
12001b : 



3D 2 N 

io- 24 /? 3 



, m 



1.3 x 10- U R 

tdyn = = , Yr 

M = 5 x 10~ 31 m p r m n e V, M e 



(1) 

(2) 
(3) 



E = 2.64 x lO^T s R 3 n e , erg 



/, = 4 x 10 n e tdy n , m 3 s 



(4) 
(5) 



We note that our distance estimate of 60 kpc to IKT 16 may 
not be entirely accurate due to the large line of sight extent of 
the SM C, which is measured to be 4-6 kpc i n this region of the 
galaxy (Subramanian & Subramaniam 2009). Errors in the dis- 
tance measurement will also impact the measured radius of the 
SNR (oc D) and its volume (oc D 3 ), and the physical properties 
derived from the above system of equations will be affected. 

Details of the best fit spectra obtained are given in Table [3] 
and the combined pn and MOS spectra for both extraction re- 
gions are shown in Fig. [3] Due to the limited quality of the 
spectral data, shock and electron temperatures are assumed to 
be the same (valid as an evolved SNR should be approaching 
electron-ion equilibration). The metal abundances in the fits are 
fixed at 0.2 solar for the same reason. Since lines, whose pres- 
ence is required to fit the ionization age, are not easily discern- 
able in the X-ray spectrum, this parameter was calculated based 
on the above system of equations. Spectral fitting was then re- 
peated with the calculated ionization age to derive new Sedov 
parameters, and this process was repeated until convergence was 
achieved. The best fit physical parameters from the Sedov model 
fit are shown in Table [4] The best fit models found require a sig- 
nificant absorption component within the SMC in addition to 
foreground Galactic absorption. We initially left this absorption 
as a free parameter for both regions. In this case, we find the 
best fit Sedov model for the SNR to have an electron tempera- 
ture of 0.71 keV with absorption inside the SMC of 4.6 x 10 21 
cm -2 (x 1 = 1100.7/1032), outside the error bounds of the ab- 
sorption found in the bright source spectrum. Setting the SNR 
absorption at the level found for the bright source does not sig- 
nificantly worsen the fit (x 2 = 1103.2/1033), but it does alter 
the best fit temperature (see Table |3j. We fix the absorption in 
the SNR to the level of the bright source, as is expected if the 
two emission regions are physically connected, and derive the 
Sedov parameters for the SNR based on this fit (see Table [4}. If 
we instead fix the absorption for the bright source at the level 
found for the SNR, we obtain a significantly worse fit for this 
component (x 2 - 164/137). In this case, we find the power-law 
index of the source to be somewhat steeper (F = 1.74 ± 0.08). 
For both components, the average flux across all observations 
is measured to be 2.7 x 10~ 13 erg cm -2 s -1 (0.5-10 keV), cor- 
responding to an unabsorbed luminosity of 1.6 x 10 35 ergs -1 . 
The flux for both components is found to be constant between 
observations to within 15%. There is thus no evidence for any 
long-term variability of the point source. 

The X-ray spectrum of the bright source within the SNR is 
typical of many astrophysical sources. The possible nature of 
this source will be discussed in SectionF 
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IKT1 6 spectrum Bright Source Spectrum 




Fig. 3. EPIC spectra from IKT 16 (left) and the bright source located near its centre (right). In both cases the solid line corresponds 
to the best-fit spectral model (see text and Table [3J, with the upper line showing the combined spectra from pn observations and the 
lower line combined spectra from MOS1 and MOS2 observations. The^ 2 residuals with respect to the best-fitting model are also 
shown. 



Table 3. Parameters of the best-fit model to SNR and bright source regions. 



Region 


SMC Absorption 


Power-law Cont. 


Sedov model 


Goodness 


Unabsorbed Lx" 




10 21 cnr 2 


Photon Index 


Shock T (keV) 


Ion. time (10'°cm 3 s) 


of Fit 


0.5-10 keV 






Normalization 


Normalization 




* 2 Mof 


(10 35 ergs" 1 ) 


SNR 


3.4(f) 


1.58(f) 


1.03±0.12 


6.1 


1103/1032 


1.6±0.4 






1.0 x l(r 5 (f) 


1.4 ± 0.3 x 10~ 4 








Bright Source 


3.4±0.6 


1.58±0.07 


1.03(f) 


6.1(f) 


150/137 


1.6±0.2 






3.4 ± 0.4 x 10~ 5 


3.6 x 10- 6 (f) 









° - Includes correction for spillover of bright source photons into SNR extraction area and vice versa, 
(f) - Parameter fixed for consistency between fit regions. 



3.2. Optical 

Using data from MCELS, we find that the [S n]/Ha ratio found 
in the SNR is ~ 0.25 (Fig. [1] top right), significantly below the 
standard "minimum" ratio typically found in radiative shocks. 
This implies that there is no evidence for radiative shocks around 
the SNR. In known SMC SNRs, this ratio is typically found to 
be higher than 0.4 dPavne et al.ll2007l) . 

Similar results are found from analysis of the central region 
of IKT 16. Using QFitsViewfl a one-dimensional spectrum of 
the diffuse region excluding all visible stars was created (Fig. |4}. 
Analysis required the identification of cosmic rays on the 2- 
dimensional spectra while the IRAF task splot allowed determi- 
nation of emission line counts based on gaussian fits. Standard 
errors were calculated based on an average rms noise of 600 
counts; propagation of error methods were used for the addition 
or division of line fluxes. Using this method, we find a [S ii]/Hq- 
ratio of 0.23 with a 2.1% error. It is quite possible that there are 
very faint nebulae interspersed between the stars that do not meet 
the criteria for shocked regions. We cannot exclude the possibil- 
ity of an HII region being present in the line of sight. 



4 Written by Thomas Ott and freely available at 
www.mpe.mpg.de/ ott/dpuser/index.html 
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Fig. 4. WiFeS spectra of the central region of IKT 16. 



3.3. Radio 

The available radio images show IKT 16 to be a very radio- 
weak SNR. There is evidence of faint diffuse radio structure 
throughout the remnant. The images are dominated by a source 
corresponding to the X-ray bright source, which appears to ex- 
tend a distance of ~ 30" towards the centre of the remnant. 
There is no evidence for a radio point source located within this 
emission. Table [2] shows the integrated flux density measure- 
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ments at frequencies from 843 MHz - 8640 MHz. The overall 
radio spectrum (Fig. |2) is well-described by a power-law with 
index a = -0.45 + 0.02. The radio spectral index of the cen- 
tral source is found to be flatter (-0.19 ± 0.02) than that of the 
SNR itself (-0.97 + 0.35). The SNR spectral index is unusu- 
ally steep, especially given that IKT 16 is most likely an older 
(evolved) SNR, due to its rather large size. Usually, a steep gradi- 
ent like this would suggest a much younger and more energetic 
SNR. However, in this case, the steepness can be contributed 
to missing short spacings at higher radio-continuum frequencies 
(4800 and 8640 MHz) and therefore missing flux. Specifically at 
8640 MHz (where the ATCA primary beam is -300") this SNR 
edges would be positioned close to the primary beam boundary 
where the flux tends to be significantly uncertain. We note that 
if we omit the two highest frequency data points from the SNR 
fit, the radio spectral index is -0.72, which is more in line with 
that found in other SNRs in the SMC (Filipo vic et al.ll2005l) . 



4. Discussion 

The only significant study of the X-ray properties of this remnant 
to date has been completed by VDH04. The XMM-Newton ob- 
servations covering IKT 16 since publication of this paper have 
significantly improved the available statistics. This allows us to 
examine the properties of the SNR more accurately, as well as 
enabling study of the nature of the bright source located close to 
its centre. 

4.1. Properties of the SNR 

In the analysis of IKT 16 conducted by VDH04, the bright 
source emission area was excluded from analysis, and the re- 
maining emission was fit with both a single-temperature non- 
equilibrium ionization (NEI) model and a Sedov model. As we 
have shown, simply excluding the source emission area may not 
produce an accurate measure of the properties of the SNR, since 
~ 25% of the flux detected through this method will be from 
the point source. As this source is significantly harder than the 
SNR as a whole, this unduly influences the results obtained. 
Additionally, the previous study estimates the radius of the SNR 
to be 100" compared with 128" here. The new estimate is based 
on mosaicing of all available observations, through which the X- 
ray shock boundary is more clearly defined in all directions. This 
changes the estimate of the X-ray emitting volume, which also 
affects the derived Sedov parameters. 

The Sedov properties obtained for both studies are shown 
in Table [4] The properties we find justify our use of the Sedov 
approximation, as the swept-up mass significantly exceeds the 
ejecta mass but radiative cooling has not occurred yet. The 
parameters we derive are significantly different from those in 
VDH04. The temperature we find is cooler, which we attribute 
to the different methods of dealing with the point source men- 
tioned in the previous paragraph. Along with our larger estimate 
of the radius of the remnant, this leads to a higher dynamical 
age, larger swept-up mass and lower initial explosion energy (as 
expected from EqJT]|4]>. In order to test whether this explains the 
difference in our results, we attempted to fit the SNR emission 
region with a Sedov model without accounting for spillover from 
the bright source. We find a good fit Ot 2 = 1099/1032) with a 
higher shock temperature of 1.30 ± 0.20 keV. This is consistent 
with the value found by VDH04 ( 1 .76+0.65 keV). VDH04 noted 
that the initial explosion energy appears high in comparison to 
typical SNRs, which they suggest is due to the remnant resulting 



Table 4. Physical properties of IKT 16 from best fit Sedov 
model. 



Parameter 


Units 


This paper 


VDH04 


Temperature 


keV 


1.03 


1.76 


Radius 


pc 


37 


29 


Electron Density 


cm" 3 


0.03 


0.05/0.04" 


Dynamical age 


yr 


14700 


7500 


Swept-up mass 


M G 


232 


124 


Explosion energy 


10 51 erg 


1.2 


3.4 



" - VDH04 derive electron and proton densities separately, assumed to 
be approximately the same here. 

from collapse of a massive star within a low-density circumstel- 
lar cavity. Our estimate of the explosion energy (~ 10 51 erg s _1 ) 
is much more in line with that derived for other SNRs in the 
SMC. 

Compared with other SNRs in the SMC, IKT 16 is notably 
X-ray faint and is detected at higher electron temperature ( ~ 1 .0 
keV vs. ~ 0.2 keV). While this is the largest confirmed SNR in 
the SMC there are several larger cousins in the Large M agellanic 
Cloud (LMC) (ICaiko et al]l2009t iKlimek et alJl2010h . The ra- 
dio emission from the SNR is also found to be faint (32 mJy at 
1400 MHz; Table [2]). The large size of the remnant, consider- 
ing that it is still in the Sedov phase of evolution, suggest that 
the density of the SNR environment is low. This may also con- 
tribute to its X-ray and radio faintness. The low density implied 
is consistent with a SNR explosion inside a wind blown bubble, 
which are typically found inside molecular clouds. Such sites are 
associated with massive star formation, which would increase 
the likelihood of this system being a c ore-collapse SNR an d the 
central source being associated with it. Mul ler et al.l (12010) have 
carried out a survey of CO emission across the northern SMC, 
finding a series of molecular clouds near the position of IKT 16. 
Their presence strongly supports this scenario. The apparent X- 
ray faintness of the remnant may be additionally enhanced by 
high absorption atte nuating the soft X - ray si gnal. Using a map 
of HI density from Stanimir ovic et al.l (1 19991) . the total column 
density of the SMC in the direction of IKT 16 is 5.7 x 10 21 crrT 2 . 
The value we derive from spectral fitting is «60% of this value, 
indicating that IKT 16 is located deep within the SMC. 

The soft X-ray image shows that the distribution of X-ray 
emission appears somewhat uneven inside the remnant, with 
more emission observed towards the northern shock boundary 
than to the south. The radio images also show more emission 
to the north of the remnant. The 3-colour X-ray image (Fig. fl] 
top left) suggests that the X-ray emission in the north may be at 
a higher temperature than in the rest of the remnant. This may 
be due to a combination of several factors, including fluctua- 
tions in the SMC absorption column, a higher concentration of 
SNR ejecta increasing the local metallicity or true variations in 
the plasma temperature of the ISM. Due to the limited X-ray 
statistics, it is not possible to carry out spatially resolved spec- 
tral analysis to investigate this further. 

4.2. Nature of the Bright source 

The X-ray spectrum of the bright source is consistent with that 
of several different astrophysical objects, including X-ray binary 
systems, pulsar wind nebulae and background active galactic nu- 
clei. Here, we attempt to discriminate between these possibili- 
ties. 

In fitting a power law model to the source, a significant 
absorption component is required in addition to foreground 
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Galactic absorption. This suggests that it is not a Galactic ob- 
ject located along the line of sight. 

It is possible, however, that the source is a background AGN. 
The average power-law phot on index o f a type I AGN is mea- 
sured to be ~ 1.9 dPounds et al.lll990l) . somewhat softer than 
what we see here. However, the absorption required in our spec- 
tral fit (3.4 ± 0.6 x 10 21 cm -2 ) is significantly lower than that of 
the SMC along our line of sight, making this unlikely. We there- 
fore attempted to fit a more complex model commonly found to 
fit AGN spectra: a power law p lus a 0. 1 keV black body com- 
ponent dTurner & Po unds 1989). Here, we find a reasonable fit 
with an extra absorption component of 3 ± 1 x 10 21 crrT 2 . The 
X-ray flux detected in this model is ~ 3 x 1CT 13 erg cirT 2 s _1 . 

Studies of type I AGN typically show the fluxes in X- 
ray (0.5-10 keV) and optical bands to lie within approxi- 
matel y an order of magnitude of each other (Maccacar o et al.l 
Il988h . For a source with this observed X-ray flux, we there- 
fore expect there to be an optical counterpart at a V mag- 
nitude of between 14.8 and 20.8. We s earched the Optical 
Gravitational Lensing Experiment (OGLE) (lUdalski et al.f l 998) 
and Magellanic Clouds Photometric Survey dZaritskv et al.l 
|2002|) to try to find the optical counterpart. The closest match 
found to the X-ray position (searching down to a V band op- 
tical magnitude of 21.5) is at a distance of 2.5" , 2 times the 
combined errors in optical and X-ray positions. However, the 
two closest optical sources in the catalogue have optical colours 
consistent with stars within the SMC. The closest source with 
optical colours atypical for a main sequence star is at a distance 
of 3.2" , still within 3 sigma of our X-ray source position. This 
is a possible optical counterpart to our source. 

We can also estimate the probability of finding a background 
AGN coincident with the S NR position. Studie s of the log N-log 
S distribution of AGN (e.g. iRosati et ai]|2002l) show that above 
our estimate of the unabsorbed X-ray flux of this object (~ 4 x 
10~ 13 erg cirT 2 s _1 , 0.5-10 keV) we expect to find 3 AGN per 
square degree of the sky. The probability of finding one such 
source within 30" of the centre of IKT 16 is 0.07%. Although we 
are unable to rule out an AGN origin for this source, we conclude 
that it is statistically unlikely. We assume from this point that the 
source is located inside the SMC. 

As shown in ^3] the SNR shell is well-described by a circle 
of radius 128" . Presuming that the well-defined profile of the 
western shock boundary can be used to find the centre of the 
remnant, the bright source is located as 30 + 8" (a transverse 
distance of 8 + 2 pc) away from this point. Providing it is asso- 
ciated with the supernova event, it has traversed this distance in 
the dynamical age of the remnant (14700 yr), at an average speed 
of 580 ± 100 km s" 1 . The direction of travel is not necessarily 
perpendicular to us, so this is a lower limit. However, unless the 
direction of travel is significantly out of this plane, the velocity of 
the object will not be much greater than this value. Asymmetry 
in SNR explosions typically gives the resulting compact object 
a rand om space velocity of 350 - 550 km s _1 dLvne & Lorimerl 
1 19941) . so it is reasonable to suggest that the bright X-ray source 
is related to this object. There are two known objects which may 
produce the emission observed: a pulsar wind nebula or a micro- 
quasar. 

Pulsar wind nebulae (PWNe) are generated by dissipation of 
a pulsar's rotational energy through relativistic winds, producing 
synchrotron emission visible in radio to X-ra y bands. For a re- 
view of the properties of these systems, see iGaensler & Slane 
d2006l) . ~ 50 PWNe have been observed in the Galaxy and 
the LMC, and the majority are assoc iated with known SNRs 
("Composite" SNR. iKaspi et al. 2006). The X-ray spectrum of 



known PWNe is well-describ ed by a power-la w of index rang- 
ing between F ~ 1.1 — 2.3 (Chen sTet al.ll2004l) . consistent with 
the power law fit found here. The measured X-ray luminosity of 
1.6 x 10 35 erg s _1 is also sim ilar to that of known PWNe in the 
Galax y and the LMC ( e.g. IGaensler et"al]|2003bt iPorquet et all 
2003; Slan e et al. I l200l . In our observations, the maximum ex- 
tension of the X-ray source was found to be 6" , which cor- 
responds to a radiu s of ~ 1 pc at the distance of the SMC. 
Gelfa nd et al. I d2009h find that in the spherically symmetric case, 
a pulsar wind nebula inside a Sedov-phase remnant at an age of 
15 kyr should be undergoing compression due to the SNR re- 
verse shock to a radius of < 1 pc. This would be detected as a 
point source in XMM-Newton observations, consistent with our 
measurements. The radius measured in the Chandra observation 
is 0.3 ± 0.2 pc, consistent with this scenario. In this case, the 
movement of the source away from the centre of the SNR shell 
implies that spherical symmetry no longer holds, but as the dis- 
tance travelled is only ~ 10% of the shell diameter, this scenario 
may still be valid. Unfortunately, the limited spatial resolution 
of XMM-Newton and the distance to the SMC (60 kpc) preclude 
identification of a central pulsar, which would confirm the iden- 
tity of this source as a pulsar wind nebula. The relatively low 
X-ray flux also prevents us from conducting timing analysis to 
search for pulsations. 

The ATCA images in FigJT] show the source to have a ra- 
dio extent of ~ 40" x 30" , elongated in the direction of the 
centre of the remnant. This extent is significantly greater than 
that observed in X-rays. This may be due to the longer cool- 
ing time for radio emission compared to that of X-ray emission. 
The radio emission may therefore trace the path of the neb- 
ula as it moves from the centre of the SNR , as is observed in 
"relic" PWN (e.g. Vela X, iFrail et al.l[l997h . The radio spectral 
index of this source is measured to be ~ -0. 2, consistent with 
the ty pical radio spectrum found in PWNe (IGaensle r & Slane 
2006). The absence of a radio point source in this emission rules 
out the presence of an established radio pulsar. However, con- 
firmed PWNe (e g. LMC SNR B0453-685 and L MC SNR N 206; 
IGaensler et all d2003al) : IWilliams et all (120051) : Crawford et al. 
(in prep)) have been found where the central pulsar is not de- 
tected in radio observations, so this does not discount a PWN 
origin. 

The other possibility we consider here is a microquasar trav- 
elling through the SMC, where the extended radio emission is 
due to optically thin synchrotron radiation in relativistic jets. 
The X-ray emission in this case would come from an X-ray 
binary (XRB) system with a power-law of index 1 .5-2.2, with 
radio jets observed at a spectral index betw een -0.8 and -0.4 
(e.g. lFendenl2006l: Hannikainen et al. I I2006T) . The X-ray power- 
law slope is consistent with this scenario, but the radio spectrum 
is somewhat flatter than expected (a ~ -0.2). If it is associ- 
ated with the SNR, and providing the western shock boundary 
can be used to estimate the original site of the supernova explo- 
sion, the velocity we infer for this object (580 ± 100 km s~') 
is very high compared with the average kick velocity for even 
Low-Mass X-ray Binary (LM XB) systems (180 ± 80 km s -1 ; 
iBrandt &" Podsiadlowski[ ll995l) . In fact, the majority of known 
microquasars are High-Mass X-ray Binary (HMXB) systems, 
for which t he corresponding typical kick v elocity is 50 km s _1 
(also from IBrandt & Pod siadlowski (1995)). It is possible that 
the object may be within the SMC but not associated with the 
SNR, in which case this velocity measurement is not applica- 
ble. The relative stability of the X-ray flux between observations 
is somewhat unexpected, as many X-ray binary systems will ex- 
hibit state changes on a timescale of 10 years. From the available 
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X-ray and radio evidence, we conclude that it is unlikely that the 
bright source is a microquasar. We note that no LMXB systems 
are currently known in the SMC, and they are very rare in the 
Galaxy. This also decreases the likelihood of this scenario. 

On the basis of current observations, it is most likely that 
the source is a PWN. However, we cannot entirely rule out a 
microquasar or AGN. A future deep Chandra observation will 
be able to determine the X-ray extent of this source accurately, 
and will thus conclusively determine its nature. 

5. Conclusions 

We have carried out analysis of the spatial and spectral properties 
of the SMC SNR IKT 16, using archive and proprietary data 
from XMM-Newton, radio data from ATCA and optical images 
from MCELS. The properties of the SNR are consistent with it 
being in the Sedov-adiabatic phase of evolution. The dynamic 
age of the remnant and initial explosion energy are typical of the 
population of SNRs in the SMC. 

We find the X-ray point source near the centre of IKT 16 
to be non-thermal, with an X-ray spectral index of ~ 1.6. This 
source is coincident with a radio source which extends a dis- 
tance of ~ 30" in the direction of the centre of the SNR. We 
conclude that the X-ray source is unlikely to be a background 
AGN, due to the small probability of spatial coincidence with 
the SNR and the absence of an obvious optical counterpart. Two 
astrophysical objects are considered which may explain the X- 
ray and radio emission observed: a pulsar wind nebula or micro- 
quasar. The flatness of the radio spectrum, the one-sided mor- 
phology of the radio emission and inferred velocity of the source 
strongly favour a PWN origin, which may be confirmed by a fu- 
ture Chandra observation. If this is the case, IKT 16 is the second 
largest such system (after the Vela SNR) found to date. 
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